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Abstract 


Describes a computer program (Fortran IV 
language) designed to analyze hydrographs. 
Using only input data of stage and time, the 
program automatically determines the begin- 
ning and end of each event, analyzes each 
event with respect to all its hydrologic charac- 
teristics (including a partitioning of base flow, 
interflow, and rapid flow), and then prints a 
detailed analysis of each event and a series of 
summaries of all events analyzed. The hydrol- 
ogist can influence the results only by chang- 
ing the time interval between hydrograph 
readings. The method has been tested on both 
simple and complex hydrographs representa- 
tive of runoff events in many regions of the 
United States. 


INTRODUCTION 


A hydrograph is a curve depicting a 
stream’s rate of flow throughout a given pe- 
riod of time. If the records are complete, the 
curve is continuous; this curve reveals to the 
trained observer the innermost secrets of a 
watershed’s behavior — nothing is concealed. 
Every change in the curve reflects the water- 
shed’s reaction to a multitude of internal and 
external forces. A steep rise in the curve could 
reflect the passage of a high intensity storm; a 
gradual descent might reflect the watershed’s 
high water storage capacity. Every change in 
the hydrograph can be readily interpreted i. 
the light of forces acting on the watershed. A 
study of all changes in the hydrograph is a 
key to our understanding of the functioning of 
a watershed. 

Streamflow is one of the end products of 
precipitation that falls on a watershed. Some 
of this precipitation falls directly into the 
stream channel and brings about an immedi- 
ate rise in the rate of streamflow. Similarly, 
another portion of the precipitation quickly 
enters well defined channels; it also contri- 
butes to the rapid rise in streamflow rate. 
Some of the rain enters the soil mantle, and 
thence, by circuitous paths, finally enters the 
stream channel. Another portion, pulled by 
the force of gravity, finds its way to the water 
table (the surface of the groundwater), and 
eventually flows into the stream channel. The 
shape of the hydrograph is determined by the 
amount and rate at which each of these com- 
ponents enters the stream channel. This shape 
is a unique expression of a watershed’s re- 
action to its environment. 


Although the ultimate purpose of hydro- 


graph analysis is an understanding of a water- 
shed’s behavior, this can be achieved only if 
we have a thorough understanding of the flow 
patterns that together make up the hydro- 
graph. For this reason, hydrologists have de- 
vised many methods to separate the com- 
ponents, despite the belief that “‘it is prac- 
tically impossible to measure the volumes of 
water following each path... .”” 

The first and most important step in 
hydrograph analysis is to partition the hydro- 
graph into its components: 

Total flow = Base flow+Interflow+ Rapid flow 

Base flow is the contribution from ground- 
water — water that moves under differences in 
head in a saturated zone. Interflow is water 
that percolates through the soil mantle and, 
after some delay, enters the stream. Rapid 
flow (also called surface flow) is water that 
enters the stream channel with a minimum 
delay; it includes channel interception. The 
sum of interflow and rapid flow are desig- 
nated as stormflow. 

The techniques now employed in hydro- 
graph analysis are highly subjective. Moreover, 
the analyst starts his analysis at the end of the 
storm hydrograph and proceeds backwards in 
time towards the beginning of the storm 
hydrograph.’ Linsley and others discuss thor- 
oughly the procedures and arbitrary rules that 
are used in hydrograph separation. 

This paper presents a computer program 
designed to analyze automatically hydrologic 


‘Ray K. Linsley, Jr., Max A. Kohler, and Joseph 
L. H. Paulhus. P. 387, in: Applied Hydrology, 689 p. 
New York: McGraw-Hill Co. 1949. 


2 Ibid, Ch. 15. 


events recorded on a hydrograph. These 
events could result from either, or a combina- 
tion of, rainstorms, snowmelt, or excessive 
evapotranspiration. The analytical technique 
is completely automated; the results are de- 
pendent entirely on the flow rates shown on 
the hydrograph. The hydrologist can influ- 
ence the results in only one way: by increas- 


ing or decreasing the number of sampling 
points along the hydrograph. 

The technique used in the computer pro- 
gram was tested on hydrographs of typical 
runoff events in many regions in the United 
States; it was considered successful because 
the computed separation curves crossed the 
recession limbs of the hydrographs at their 
points of inflection. 


BASIC ASSUMPTIONS 


Two basic assumptions underlie the com- 
putations: that both increment in base flow 
and the ratio of surface runoff to interflow 
depend essentially on the hydrograph’s rate of 
rise. A rapid rise indicates a proportionately 
smaller increment in base flow than does a 
moderate rise. On the other hand, the ratio of 
surface runoff to interflow is large when the 
stream rises rapidly; this ratio is small when 
the stream rises slowly. 


We also assume that (a) surface runoff 
(including channel interception) ends at the 
first inflection point following peak flow, and 
that (b) interflow ends at the lower inflection 
point on the recession limb of a simple hydro- 
graph. Our research was directed towards de- 
fining equations that — if based on the first 
two assumptions — would meet the criteria of 
the third and that would be applicable to any 
hydrograph. 


ESSENTIAL ELEMENTS 
OF THE SEPARATION TECHNIQUE 


The heart of this program consists of a se- 
ries of calculations to determine the consecu- 
tive positions of points along two curves: the 
base flow curve (which separates base flow 
from stormflow) and the interflow curve 
(which separates interflow from rapid flow.) 
Figure 1 is used to illustrate the parameters 
entering the calculation of points along both 
the base flow and interflow lines. 


Base Flow Computations 


The rate of flow at any point (P;) on the 
base flow line is equal to the sum of the rate 


at a prior time (es ) plus an incremental val- 
ue (Uj): 

et SOG (1) 
where 
U; = exp[-(1.0+sin A)]x(At x AR x tan B)” (2) 


where (see fig. 1) 
AR=R;-Pj., 
pea 
A = arctan (AR/At) 
B=(H xA -AR)/H 
H =(AR* + At? )”% 


Oe 


Figure 1. — Parameters used to locate points 
along the base flow line. 


Each successive point on the hydrograph is 
the basis for a new calculation of P;. The 
event ends when the value of a point on the 
base flow line equals a rate value on the 
hydrograph trace; i.e., Pj=R;. The next 
hydrograph event begins when the hydro- 
graph starts to rise again. At that point, 
Ri, = E. ; 


Interflow Computations 


The curve that separates interflow from 
surface flow (which includes channel intercep- 
tion) is similarly computed from point to 
point. Any point S; on the interflow curve is 
equal to the value at a prior time (Sj, ) plus 
an incremental value (W)). 

Thus, S;=S;.,+ Wi (3) 
where 


W; =[(At? + AS?)? x AS]/(At+KxAS) (4) 


where 
At=;- tie, 
AS = R; - Si-, 
K = 3.0 


Surface runoff ends when the curve inter- 
sects the hydrograph; i.e., S; = Rj. 


DISCUSSION OF TECHNIQUE 


Table 1 shows the results of analyses of 
hydrographs of representative storms in sever- 
al regions of the United States. The analyzed 
hydrographs include both simple and complex 
hydrographs. Source for the basic hydrograph 
data is Hobbs and Burford.’ In all cases, the 
analyses were performed on data expressed in 
cubic feet per second (c.f.s.). Data originally 
expressed in inches per hour were converted 
to c.f.s. Figures 2 and 3 show the separation, 
respectively, of a simple and complex hydro- 
graph into its components. The figures show 
the total flow hydrograph, and the computed 
interflow and base flow rates. 

Normally, on large watersheds, the time in- 
terval between sampled points on the hydro- 
graphs is one or more hours. Therefore, the 
analytical results will be independent of the 


3 Harold W. Hobbs and James B. Burford. Hydro- 
logic data for experimental agricultural watersheds in 
the United States, 1963. USDA Agr. Res. Serv., Misc. 
Pub. 1164. 1970. 


hydrologist. On very small watersheds, the 
hydrologist tries to reproduce the shape of 
the hydrograph by sampling at small and vari- 
able increments of time because the smaller 
the time interval used the more accurately 
can he reproduce the hydrograph’s original 
shape. Thus, on very small watersheds the 
hydrologist can influence the results by in- 
creasing or decreasing the number of sampling 
points used to reproduce the shape of the 
hydrograph. 

Our analyses show that when we reduced 
the number of sampling points 75 percent, we 
changed the computed results by significant 
amounts. The maximum change was observed 
on a small watershed (0.07 km.” ) from which 
all runoff ceased after 1.75 hours: computed 
total flow was increased 0.457 cm. (28 per- 
cent); stormflow increased 0.508 cm. (28 per- 
cent); and stormflow runoff duration was in- 
creased 6.6 minutes. Consequently, the 
hydrologist who analyzes the hydrographs of 


FLOW RATE (Inch/hour) 


Table 1. — Analysis of hydrographs of streams in different regions of the United States’ 


Vero Beach, Blacksburg, N. Danville, Hastings, Colby, Waco, Reynolds, 
Florida Virginia Vermont Nebraska Wisconsin Texas Idaho 
Area (km.? ) 255.5 6.0 Talal 8.3 1.4 223 233.4 


PERCENT OF TOTAL FLOW 


Base flow 32.47 14.80 53.16 9.07 22.74 5.52 13.36 
Stormflow 67.53 85.20 46.84 90.93 77.26 94.48 86.64 


PERCENT OF STORMFLOW 


Interflow 65.38 89.99 90.62 66.58 90.41 89.63 93.49 
Surface flow 34.62 10.01 9.38 33.42 9.59 10.37 6.51 
Before peak 44.19 32.47 37.81 39.92 14.59 26.71 14.68 
After peak 55.81 67.53 62.19 60.08 85.41 72.39 85.32 


DURATION (HOURS) OF STORMFLOW 


Time to peak 48.00 1.00 10.08 Sulai Arit2, 1.23 18.00 
Recession time 56.64 3.36 24.88 5.91 8.87 3.89 38.60 
Total time 104.64 4.36 34.96 9.08 12.99 toys 56.60 


' Source of data: ibid, text footnote 3. 
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Figure 2.— Flow components in a simple hydrograph recorded at Riesel, Texas. 
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Figure 3.— Flow components of a complex hydrograph recorded at Hastings, Nebraska. 


very small watersheds should make every ef- 
fort to select points that reproduce as ac- 
curately as possible the true shape of the 
hydrograph. 

Equation (2) includes two angle measure- 
ments (A and B). Accordingly, it is apparent 
that the results of any analysis will depend on 
the units of measurement, as shown in table 2. 

Normally, the analysis would be performed 
on data expressed in a unit suitable for the 
watershed: c.f.s. for low water-yielding water- 
sheds, and inch/hour for high water- 
yielding drainages. Certainly, a great deal of 
accuracy is lost when calculations are per- 
formed on angles measured in minute frac- 
tions of a radian. For example: suppose 
streamflow from a 580-acre (2.346-km.’ ) 
watershed is expressed in inches. If the in- 
crease is 0.0005 inch in 0.1 hour, angle A is 
approximately 0.005 radian; if expressed in 
c.f.s., the angle is 1.24 radians. The effect of 
loss in accuracy as well as the unit of measure- 
ment is fully reflected in the last column of 
table 2, which shows figures for surface run- 
off based on data expressed in inch/hour. 

Equation (4) is used to partition stormflow 
into surface runoff (including channel inter- 
ception) and interflow. The equation contains 


a constant (K) to which we assigned a value 
of 3.0. The exact value of K can be deter- 
mined only if precise data from controlled 
laboratory or field experiments are available, 
but we do not have such data. 

A change in the value of K results in sig- 
nificant changes in the analytical results (table 
3). Increasing the value of K yields (a) larger 
computed values for surface runoff; (b) a cor- 
responding decrease in interflow; and (c)a 
larger computed value for duration of surface 
runoff. Our analyses for large forested water- 
sheds show that the computed period for sur- 
face runoff duration is extended by 1-3 hours 
if the value of K is increased from 3.0 to 4.0. 

We believe that the value of K is related to 
the proportion of the watershed represented 
by the main channel and its tributaries. For 
any particular watershed, the value of K 
might well change with either or both the sea- 
son and the characteristics of the storm. This 
concept of K agrees with the concept of a 
variable source area.* 


4 John D. Hewlett and Alden R. Hibbert. Factors 
affecting the response of small watersheds to precipi- 
tation in humid areas. P. 275-290, in: Forest Hy- 
drology, edited by William E. Sopper and Howard W. 
Lull. New York: Pergamon Press, Inc. 1967. 


The assigned value of K = 3.0 is acceptable the precise partitioning of stormflow into sur- 
if the hydrologist is mainly interested in the face runoff and interflow will depend on the 
comparative behavior of several watersheds or exact value of Kas determined under controlled 
in the reaction of a single watershed to storms experiments. 
of different magnitude. On the other hand, 


Table 2. — Effect of unit of measurement on results of analysis 


Basic hydrograph unit Basic hydrograph unit 


Location Gs! inch/hr. c.f.s. inch/hr. 
Stormflow as % Surface runoff as % 
of total runoff of stormflow 
Vero Beach, Fla. 67.53 69.80 34.62 0 
Blacksburg, Va. 85.20 60.01 10.01 7.89 
N. Danville, Vt. 46.84 82.54 9.38 0 
Hastings, Neb. 90.93 87.19 33.42 8.14 
Colby, Wisc. TU26 75.12 9.59 0.78 
Waco, Tex. 94.48 66.30 10.37 ITS) 
Reynolds, Ida. 86.64 59.31 6.51 0.37 


Table 3. — Effect of K (Equation 4) on surface runoff as percent of storm- 
flow during October-November 1968, in the East Fork Horse 
Creek watershed, Idaho 


Event No. Stormflow K=2 K=3 kK=4 K=65 
Cubic ft./100 
19 6.7 Bid 8.60 11.33 13.96 
4 TRAD) 170 DP ASYS) 3.40 4,24 
3 68.8 DWN 3.34 4.55 5.81 
23 98.9 3.64 5.36 7.10 8.90 
ak 202.6 0.26 0.38 0.49 0.60 
alts) 308.4 1.54 2.44 3.44 4.53 
34 613.2 11.93 18.53 24.63 SONS 
28 1873.1 a9 19.65 27.05 33.63 
16 3898.0 11.89 19.95 27.10 33.24 
6 


—_ 
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THE COMPUTER PROGRAM 


The program consists of three sections: 
Input of basic data (lines 1 through 12); anal- 
ysis of data (lines 13 through 129); and output 
of results (lines 130 through 183). 


Input of Basic Data 


The program is designed to accommodate 
data for 1 year, providing the time interval be- 
tween readings is not less than 1 hour and the 
number of storm events is not more than 300. 
The dimension statements should be changed 
if the number of data points exceeds 9000 or 
the number of anticipated storm events ex- 
ceeds 300. If either of these conditions ap- 
plies, the author suggests that less than a 
year’s data be analyzed at one time. 

The “read’’ statement should fit the for- 
mat in which the input data are punched on 
the cards. In this program, there is no “‘read”’ 
statement for hour, because all hydrograph 
readings were taken at 1-hour intervals. If a 
variable time interval characterizes the input 
data, then a new “read”’ statement should be 
made. 


Analysis of Data 


1. Volume parameters 
Total flow 
Base flow 
Storm flow 
Interflow 
Rapid flow 
2. Time parameters 
Hour at which event starts 
Hour at which event peaks 
Hour at which event ends 


Number of hours from start to peak 

Number of hours from peak to end 
(recession time) 

3. Rate parameters 

Initial streamflow rate 

Peak rate 

Final rate 
4. Angle parameters (in radians) 

Rising angle (formed by intersection of 
horizontal line and line connecting ini- 
tial and peak rates) 

Recession angle (formed by intersection of 
a vertical line and line connecting peak 
and final rates) 

Top Angle (formed by intersection of line 
connecting initial and peak rates with 
line connecting peak and final rates) 

Base angle (formed by intersection of hori- 
zontal line and line connecting initial 
and final rates) 

5. Line parameters 

Base line (length of line connecting initial 
and final points of the storm hydro- 
graph) 

Perimeter (sum of lines connecting initial, 
peak, and final points on the storm 
hydrograph.) 

6. Ratio parameter 

Ratio of perimeter to natural logarithm of 
1.0 + stormflow) 

The relation between stage and flow rate 
(cubic feet per second) is unique to each type 
of waterflow-measuring installation. Hence, 
the user of this program must supply the cali- 
bration equation that applies to the installa- 
tion. This is done in lines 23 and 24, and also 
in lines 38, 39, and 56. In addition, the user 
must indicate the area of the watershed in 
acres (line 164). 


Output of Results 


Summaries of results of all analyzed storms 
are tabulated in a series of tables (see Ap- 
pendix 1) and also on cards. (See Appendix 2, 
lines 173 and 174, for format.) 

Appendix tables 1 through 5 are examples 
of the program printout. Tables 2-5 are sum- 
mary tables of all the storms analyzed. Table 
1 (Appendix 1) is a detailed analysis of each 
individual event, and requires some explana- 
tion. The sample printout is an analysis of a 
short storm during October 1968 on the East 
Fork of Horse Creek in north central Idaho. 

The storm started at hour 347 when 
streamflow rate was 4.139818 c.f.s. (Six deci- 
mals are printed to enable the hydrologist to 
analyze any portion of the event without loss 
of accuracy due to rounding-off procedures.) 
The storm ended at hour 361.92 when 
streamflow rate was 4.335419 c.f.s. Column 2 
shows the rate of flow (c.f.s.) at hourly inter- 
vals along the base flow line — the line sepa- 


rating base flow from stormflow. The third 
column (“runoff rate’’) shows actual stream- 
flow rate (c.f.s.) at hourly intervals. This 
figure is the input data (stage-in feet) con- 
verted to c.f.s. The last column (interflow 
rate) shows the rate of flow (c.f.s.) at hourly 
intervals along the line separating interflow 
from rapid flow. As long as values in this col- 
umn are smaller than those in the “runoff 
rate’? column, the rapid flow separation curve 
has not crossed the hydrograph. When rapid 
flow has ended, then the figures in the last 
column are identical with those in the “runoff 
rate’? column. In this example, rapid flow 
ended between hours 353 and 354. Total flow 
from hour 847 to hour 361.92 was 
68.413391 c.f.s.-hours (column 5), of which 
63.696930 was base flow (column 4) and 
4.716461 was stormflow (column 7). Storm- 
flow amounted to 6.89 percent of the total 
flow (column 6). Other figures in table 1 as 
well as in the summary tables are self ex- 
planatory. 


USES FOR THE COMPUTER PROGRAM 


Hydrograph analysis is a difficult and 
time-consuming operation. The inherent dif- 
ficulties have forced hydrologists to limit 
their analysis to a few extraordinary hydro- 
logic events, or to events that they believe 
represent average conditions. 

By using this computer program, the 
hydrologist can analyze all streamflow 
events — rainstorms, snowmelt events, and 


even daily events resulting from excessive 
evapotranspiration. Any event, regardless of 
its cause, can be analyzed automatically, and 
can be recorded in a condensed form. The 
computer program enables the hydrologist to 
readily compare the numerous elements of 
hydrologic events occurring on two or more 
watersheds, because they are recorded in 
orderly form. 


APPENDIX 1 


Tables 1 through 5 are examples of the 
printed output of the program. Table 1 is an 
analysis of storm No. 6. Tables 2 through 5 are 
summary tables of 34 storms, which include 
storm No. 6. 
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APPENDIX 2 


The Hydrograph Analysis Program 
Written in Fortran IV language 
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0001 


0002 
0003 
0004 
0005 


0006 
0007 
0008 


0009 


0010 


0011 
0012 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0C25 
0026 


OC27 
0028 
0029 
0030 
C031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
00467 
0048 
0049 
0050 
0051 
0052 
0C53 


GGG). 


DIMENSION S( 9000) sHOQUR( 9000) ,QFLL( 300) »TS(300),TP( 300), TE(300), 
22S (300) ,RE( 300) 0R1 (300) ,DR2(300),0T1 (300) »0T2(300), RP (300), 
3UP (300) ,DN( 300) »PERIM(300) »A(300),PERFL(300) »A2(300),ARISE(300)y 
4AGLDN( 300) »BASE ( 300) 
DIMENSTON QFLPK(3C0) »QFPPK(3C9) 
DIMENSTON TOTTM( 300) ,QINCH(300) 
DIMENSION SUBFF (30C) »,RFLL(300),SLBIN(300) »RINCH( 300) 
DIMENSION PTFLO( 300) »PRFLO(300) »PQBPK( 309) 4PCAPK (300) 
THF PROGRAM ACCOMODATES 300 STORMS. IF A LARGER NUMBER OF STORMS IS 
ANTICIPATED, THE DIMENSION STATEMENT ABQVE SHOULD BE CHANGED, OR 
THE NUMBER OF DATA 8EING ANALYZED AT ONE TIME SHOULD BE REDUCED. 
DO 7 J=1,9000 
S(J)=-1.0 
ENDS=-1.0 
IN THIS PROGRAM,S=STAGE.THE CONVERSION FROM STAGE TO CFS TAKES PLACE 
BELCW»UNDER 'O= ', AND 'RS(K)= *.CONVERSTON CONSTANTS SHOULD BE CHANGED 
TO FIT THE WEIR BEING USED. 
READ (5y1,ENCD=8)S 
THE FOLLOWING FORMAT STATEMENT MUST FIT THE FORMAT IN WHICH THE INPUT 
DATA FOR STAGE ARE PUNCHED. 
FORMAT (10X,12F4.34922X) 
IF HYDROGRAPH READINGS ARE TAKEN AT INTERVALS OTHER THAN ONE HOUR, 
THEN A SEPARATE "READ* STATEMENT SHOULD BE MADE. 
NO 2 J=1,900C 
HOUR(J)=J 
J=1 
K=1 
START=S( J) 
IF(S(J)eLE-ENDS) GO TO 8001 
IF(S(JS+1).GT/START) GO TO 100 
LECSCSP INSEE START) S=Jer 
GO TO 4 


100 WRITE( 6.5) K 


FORMAT(*1*',50X,*TABLE L*/15X,'STCRM NUMBER', 14) 

TS(K )=HOUR( J) 

Q=1.522*(8.0**0.026) 

RS(KI=S( J) **0*32.0 

WRITE (656) TS(K),RS(K) 

FORMAT(2X%,"STORM BEGINS AT'sF10.293X,*INITIALRUNOFF RATE IS*,F12.6 
2/2X_,'TIME BASFFLORATE RUNOFRATE BASEFLO TCTALFLO STORMFLOZ 
3 STORMFLO INTERFLO RAPIDFLO INTERFLORATE!) 

U=0 
RC=0 

UC UM1L=RS(K) 

DELAF=0 
QFLO=0 

RP(K)=RS(K) 

TP(K)=HOUR( J) 

U2=0 

SUBFL=0 
RFLO=0 

UC UM 3=RS(K) 

R1=S(J)**0*32.0 

R2=S(J+1) **0*32.0 


200 OITFR=ABS(P 2-UCUM1L) 


DIFT=HOUR( J+ 1I-HCUR( J) 
T=HOUR (S41) 
H=SQRT(NIFR*e*24D TFT #2) 
IF (CH eLToNIFR)HEDIFR 
[F(R LT NIFTIHSOIFT 

RC a(R24R1)/2.0®0] FT+RC 
ALPHA®ATAN(DIFR/DIFT) 
RETAs(H*ALPHA=DIFR)/H 
IF(RETALF.O) BETA=C 
CSA21.0+SIN( ALPHA) 
tEXP(=C SA)Y*SORT(OTFR#NIFT*# TAN (BETA) ) 
UC UM 2=U#UC UML 
IF(R2.GToPP(K))TP(K)=T 
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0054 IF(R2.GTRPIK)) RP(K)=R2 

0055 IF(S(J#2).LE-ENDS) GO TO 8001 

0056 R3=S(J+2)**QO*32.0 

0057 DIFT2=HOUR( J#1)-HOUR( J) 

0058 YIFR 2=R 2—-UCUM3 

ocs9 IF(DIFR2.LE.9) GAMMA=0.0 

0cé60 IF(OIFR2.LE.0) GO Te 15 

0061 H2=SQRT(OIFR2**2+01F T¥*2) 

0062 GAMMA=H2/(DIFT2+3.0*DIFR2) 

0063 15 U2=GAMMA*D [FR2 

0064 UC UM 4=UC UM 3+U2 

océ5 IF (UCUM4.LE.’UCUM2) UCUM4=UCUM2 
0066 IF (UCUM4.GE.R2) UCUM4=R2 

0067 IF (UCUM2-R 2) 10,11,12 

0068 1C DELAF=DELAF+4( UC UML*#UCUM2)/2.0*DIFT 
0069 PERC=100.9-DELAF*¥1C0.0/RC 

0070 QFLO=RC-DELAF 

9071 SUBFL=SUBFL#( UC UM34UCUM4-UCUM1-UCUM2 )/2.0*DIFT2 
0072 RFLO=QFL9- SUBFL 

0073 TECTCEO. TP(K)) QPEPK (UK) =CELE 

OC74 WRETLE C6513) Ts UCUM2 »,R2,DELAF yRCyPERC,QFLO,SUBFL»,RFLO,UCUM4S 
0075 Ns) FORMAT(1X,F7.2,9F11.6) 

0076 R1=R2 

OoOo77 R2=R3 

0078 UC UM 1=UC UM2 

0079 UC UM3=UC UM4 

0080 J=JF1 

0081 IF(S(J).LE.-ENDS) GO TO 8001 

COR2 GC TO 200 

00F3 11 CELAF=DELAF+(UCUM14UCUM2) /2.0*DIFT 
0984 PERC=100-0-DELAF ¥*10C.0/RC 

0085 QFLO=RC-DELAF 

0086 SUBFL=SUBFL4(UC UM3+UC UM4-UCUM1-UCUM2 )/2.0*DIFT2 
0087 SUBFF(K)=SUBFL 

o0geR RFLP =QFLO-SUBFL 

0089 RFLL(K)=RFLO 

0090 QFPPK(K)=QFLL(K)-OFLPK(K) 

0091 WRITE (6413) Ty UCUM2,R2,DELAF,RC,PERC,QFLOeSUBFL»,RFELOyUCUM4 
0092 RE(K)=R2 

0093 TE(K)=T 

0094 QFLL(K)=QFLO 

99095 PIFLO(K)=SUBFF(K)*100.-0/QFLL(K) 
9096 PRFLO(K)=100.0/QFLL(K) *RFLL(K) 
0097 POBPK(K)=1C0.0/QFLL(K) #QFLPK(K) 
00°8 PQAPK(KI=1C0.0/QFLL(K) *QFPPK(K) 
0099 K=K+l1 

0100 J=J+1 

0101 PECS CI) be. ENDS) GG TC 8001 

0102 GO TO 4 

0103 12. RB1l=(R2-R1L)I/DIFT 

9104 B2=U/DIFT 

0105 TIME=(-R14+UCUM1) /(B1-B2) 

01906 TIME X=HOUR( J) +TIME 

9107 UC UM2=UC UM14+B2*TI ME 

0108 DELAF=DELAF+4( UC UM1#UCUM2)/2.0#TIME 
0109 XR =UC UM2 

0110 RC=RC-(R24+XR)/2.C*(DIFT-TIME) 

9111 PERC=1LOC.C-DELAF*1CO.O/RC 

0112 QFLO=RC—-DEL AF 

9113 OQFLL(K)=OFLO 

0114 OFPPK(K)=QFLL(K)—-QFLPK(K) 

9115 SURBFL=SUBFL4¢(UC UM3-UCUML)/2.C¥#TIME 
0116 SUBFF(K)=SUBFL 

0117 PFLN=QFLO-SUBFL 

0118 RFLL(K)=RFLO 

0119 WRITE (6,13) TIMEX, LCUM2,XR»CELAF,RCyPERC,OFLC,SUBFL» RFLO,UCUM2 
G120 RE(K)=XR 

0121 TE(K)=TI MEX 
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0122 
0123 
0124 
0125 
9126 
0127 
0128 
0129 
0130 
0131 
0132 
0133 


0134 
0135 
0136 


0137 
0138 
0139 
0140 


0141 
0142 
0143 
0144 
9145 
0146 
0147 
0148 
0149 
0150 
0151 
0152 
0153 
0154 
0155 
0156 


0157 
0158 
0159 
0160 
0161 
0162 


0163 


0164 
0165 
0166 
0167 
0168 
0169 


0170 
0171 
0172 


QiGuG 


PIFLO(K)=SUBFF(K)*1CO.C/GFLLI(K) 
PRFLO(K)=1C00.0/QFLL(K) ¥RFLL IK) 
PQBPK(K)=100-9/0FLL(K) *QFLPK(K) 
PQAPK(K)=100.0/QFLL(K) *OFOPK(K) 
J=J+1 
K=K+1 
IF(S(JS).LEENDS) GO TO 8001 
GC TO 4 

8001 WRITE( 6, 8110) 

8110 FORMAT('1",50X,*TABLE 2°) 
WRITE( 6, 80C3) 


8003 FORMAT(20X_, "SUMMARY OF STCRM DATA*/1X-* AG TIME START TIME PE 
2AK TIMEEND R START R PEAK R END STORMFLO INTERFLO RA 
3PIDFLO") 

II =K-1 
CO 8100 J=1,II 
WRITE(6, 8908) Je TSOI) o TPIS) pe TEC I)» RSC IV RPC I) eREC IID, OFLL ISD, 


2SUBFE( JV eP FERC 
8008 FORMAT(1Xs41359F 10.4) 
B190 CONTINUE 

WPITE( 6, 8902) 

8CD2 FORMAT (1 50Xe* TABLES */1X_9 "STORM e2Xy*STORMELG!® 9 29X_y* TIME’, LEX, 'P 
Z2ER IMETER/'44X_,"RECESSICON's * PASE LINE' 3Xy"RISING', 3X, "FALLING" /1X 
34°NO",3X,"BASE LINE", 3X9 '"PERIMETFR' 4X "TOP ANGLE" ,3X_,'TO PEAK"g 5X 
4, "STORMFLO%,2X_ "LN(L+#STORMFLO) " y4X5"TIME' sy 7X p* ANGLE" » 3X9 *LIMB ANGL 
SE UX SETA SANGEE ©) 

0C e200 J=1.II 
OR 1 JV=RPCJI-RS( J) 
M220 J)=RP( JI—RE (J) 
PT LM = TEE) Sie) 
NT 2 IIHT CS =— TP) 
UP ( JV=SORT(DR1( J) **24+NT1 (J) * #2) 
DNOSV=SORT(OR2A( J) **240T2(5) #42) 
PASEC( JV=SORTI( (FEU ISY-“RSCI)) 824 0TE (I) -TSC Ue 82) 
PFRIM(J)=UP( JDFENN( JS) + BASF (OS) 
A(JIV=ATAN(DTICJI/IRI(C SDD FATAN(NT2(JSI/OR2059) 
A2(JISATAN((RE( SIKRSOC ID) ACTF CJ) HTS OSD) 
ARTSE( JS) =ATAN(ORI( JI /0T1L¢08)) 
AGLDN(J)=ATAN(OT2(J)/0R2(45)) 
PERFL(J)=PFPIM( J) /ALOG((CFLL(J)41.0)) 
TOTTM( J) =TE( JI-TS( J) 
WRITE(6,8009) J+BASE(S) sPERIM(J) »A( JS) eOT1( JS), QFLL(J),PERFL(S Dy 
2 DT2(J5),A205),ARI SE( JS) »,AGLON(J) 

8009 FORMAT(1Xs13592X,8F12e6,2F926) 

8200 CONTINUE 
WRITE( 6, 8120) 

6120 FORMAT('1',50X,*TABLE 4°) 
WRITE (6,8703) 

8703 FORMAT( 1X,"STCRM NO STCR¥FLC INTERFLO RAP IDFLO QINC 
2H INTERFLINCH RAPINCH TOTAL TIME QFLO TO PEAK QFLOAFTER P 
ZEAK*) 

AINCH IS A CONSTANT FOR CONVERTING CFS INTO INCHES 
AINCH=3609.0*6490.0¥*12.0/(5280.0**2) 
ACREK IS THE AREA OF THE WATERSHED BEING STUDIED IN ACRES 
THE VALUE OF "ACREK* MUST BE CHANGED TO FIT THE WATERSHED 
BEING STUDIED. 
ACREK=3561.0 
NC 8720 J=1,IT 
QINCH( J) =QFLL(JS) *¥AINCH/ACREK 
SUBIN( J) =SURFF( J) ¥AINCH/ACREK 
RINCH(J)=RFLL(JS) *ATNCH/ACREK 
WRITE (658700) JeQFLL(J) »SUBFF( J) »RFLL( J), QINCH{ JS) ,SUBIN(J Dy 
2RINCH( J) ¢TOTTM( JS) »QFLPK( JS) »CFPPK(J) 

8700 FORMAT(1X,134295X,9F 12.6) 

8720 CONTINUF 
07 8600 J=1,II 
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0173 


0174 
QET5 
0176 
OLTT 
0178 


0179 
e120 
9181 
0182 
0183 


8600 WRITE (7,8000) A2( J) ,sARTSE (J) -A0 JS) sAGLDON( J), 0T1L(5),0T2(5),TOTTM (JS) 
29OR LIS) eDR20 I) ,QFLLO IS) » PERIM( J) pPERFLO JD eOINCH( JS) TS( IS), TEC SD, 
BRP (JS) eQFLPK( J) QFOPK (J) »SUBFF (J) sRFLL( J) 
BOCC FORMAT(5F15.6/5F15.6/3F15.692FT. 2/5 Fl2 eb el TXy"401") 
WRITE (6,813C) 
8130 FORMAT CULL 50Xs * TABLE 5°) 
WRITE (6, 7C01) 
7001 FORMAT(10Xe"SUMMARY OF STCRM FLCW CATA*/1Xe*STORM NO*,2X_,*STORM 
2FLO SINTERFLO FRAPICFLE % BEFGRE PK % AFTER PK*) 
BO SCO2 5 I= lett 
7902 WRITE(64 7003) JeQFLL(I J) eo PIFLO (J) sPRFLOC J), PQBPK( J) ,PQAPK (J) 9 J 
7993 FORMAT( LX, 1395XeF 1 32-654F13.2,14) 
9999) SToP 
END 
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Headquarters for the Intermountain Forest and 
Range Experiment Station are in Ogden, Utah. 
Field Research Work Units are maintained in: 


Boise, Idaho 

Bozeman, Montana (in cooperation with 
Montana State University) 

Logan, Utah (in cooperation with Utah 
State University) 

Missoula, Montana (in cooperation with 
University of Montana) 

Moscow, Idaho (in cooperation with the 
University of Idaho) 

Provo, Utah(in cooperation with Brigham 
Young University) 


